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Abstract 

High-frequency solutions of one or several Schrodinger-type equations are well known 
to differ very little from the plane wave solutions exp[ibi/cx]. That is, the potential terms 
impact the envelope of a high-frequency plane wave by only a small amount. However, when 
such equations are solved by a finite-difference method, the highest-frequency solutions may, 
under certain conditions, turn out to be localized. In this letter we explain this numerical 
artifact. 
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Coupled Schrodinger-type equations arise in many areas of physics, e.g., in quantum 
mechanics [U O [3l [4] and stabihty of nonhnear waves [5j. Often, the problem is posed as 
an eigenvalue problem whereby localized eigenfunctions and their eigenvalues are sought. 
The eigenfunction localization occurs due to the presence of potential-like terms. A simple 
approach of solving such an eigenvalue problem is to discretize the equations by a finite- 
difference method and then solve the resulting matrix eigenvalue problem by a commercial 
software. Then, by inspection or otherwise, one selects the localized eigenfunctions and 
their eigenvalues out of a set of eigensolutions produced by the software. 

High-frequency solutions are, typically, not sought numerically because their approxi- 
mate analytical form can be found by perturbation methods, e.g., by the Born or Wentzel- 
Kramers-Brillouin approximations. For example, for a single Schrodinger equation 

-d'^%b/dx'^ + V{x)tP = XiP, A>0, (1) 

where A ^ 1 and V{x) varies on the scale of order one and also max|y(a;)| = 0(1), the 
latter approximation yields: 

,\-V{x) 

Moreover, to accurately resolve a solution with a given frequency, one needs about 10 grid 
points per wavelength. The highest-frequency mode resolved on a grid with a step size h 
has a wavelength 2h; that is, such a mode has only 2 grid points per wavelength. Thus, 
one does not expect that such a mode can be resolved with much quantitative accuracy. 
However, one does expect that it should qualitatively look like solution its envelope 
is to be a finite constant away from the localized potential V{x) and is to have a small 
"wiggle" around the potential. 

The researcher may want to inspect the numerically obtained highest-frequency solutions 
of the eigenproblem in order to verify that at least in the high-frequency limit, his/her 
finite-difference code produces reasonable results, as described above. If the code produces 
qualitatively different profiles of high-frequency modes, the researcher may question the 
correctness of the code and search for a mistake. This can be a time-consuming task when 
a system of several coupled equations is considered. Thus, it is valuable to know what the 
highest-frequency modes obtained by a finite-difference method can look like. 

Below we show that such modes look not at all as described at the end of the para- 
graph containing Eq. ([2]). Rather, their envelopes are the Zowesi- frequency eigenfunctions 
of the potential —V{x). In short, this occurs because the finite-difference approximation to 
cP/dx'^ in ([T| evaluated on the highest-frequency carrier exjp[ikx] with k = n/h becomes 
"const — note the change of sign in front of the second derivative. Observing 

localized envelopes of the numerically obtained highest-frequency modes may be even more 
counterintuitive given that they occur for a repulsive rather than attractive potential. 
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While we noticed this fact when numerically solving several coupled Schrodinger-type 
equations, below we chose to present its explanation for a single equation ([T]), so that 
the complexity of the problem would not obfuscate the essence of the explanation. As a 
finite-difference approximation of cPip/dx'^ we use the simple central difference: 

(f^{Xn)/dx^ [tp{Xn+l) -2llj{Xn) +ll){Xn-l)) /h? , (3) 

where Xn±i = Xn^h. Using a more accurate Numerov's discretization [2] leads to the same 
qualitative conclusions. 

The following Matlab code computes the four highest-frequency modes and the corre- 
sponding eigenvalues of Eq. ([1]) with V{x) = 3 sech(0.5x) and periodic boundary conditions: 

h=0.1; x=-16:h: 16-h; N=length(x) ; 

M=spdiags(repmat( [-1 2 -1],N,1),[-1 1] ,N,N)/h-2 + diag(3*sech(0 . 5*x) ) ; 
M(l,end)=-l/h~2; M(end, l)=-l/h"2; 

[Evecs, Evals] =eigs (M,4, ' Im' ) ; absEvecs=abs (Evecs) ; 

k=l ; plot (x, sech(0 . 5*x) , ' — ' ,x, absEvecs( : ,k) /max(absEvecs ( : ,k) ) ) ; 

The envelopes of the first and fourth such modes are shown in Figs. l(a,b); the carrier is 
shown in Fig. 1(c). As we have announced above, these envelopes are localized, in contrast to 
the slightly perturbed plane waves ^ that could have been expected naively. A calculation 
that explains Fig. 1 is as follows. 
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Figure 1: (a,b) Solid lines show envelopes of the first and fourth highest-frequency modes 
obtained by the finite-difference approximation of ([I]) with V = 3sech(0.5a;). Dashed line 
show the potential The corresponding eigenvalues are shown in each panel, (c) Part 

of panel (a), magnified and displaying the highest-frequency carrier. 

As Fig. 1(c) illustrates, one can take ip{xn) = exp[iTrn] (j){xn)- Here the factor exp[i(7r//i)x„] = 
exp[i7rn] accounts for the highest-frequency carrier, while (/>(x„) is assumed to vary on the 
x-scale of order one. Substituting this ansatz into ([T|) where the second derivative is ap- 
proximated by Q, one finds: 

^<Pix) + + Vix)^ix) + Oih^) = A0. (4) 

Here x = Xn, and we have used the Taylor expansion (f){xn±i) = (t){x)^h<f)' {x)+{h? /2)(j)" {x)+ 
OQv') for the smooth envelope </>(x). Neglecting the 0(/i^)-term in Q, one sees that that 
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equation becomes 

0" + {V - AA)0 = 0, AA = A - (4//i^). (5) 

Based on ([5]), one can make the following conclusions about the appearance of the 
envelopes of the high-frequency modes of ([T]) obtained by a finite-difference approximation. 
When the potential in the original Eq. ([T]) is repulsive, V{x) > (or, more generally, 
J^^V{x) dx > 0), the envelope of the mode with the highest eigenvalue A is the ground 
state of ([5]). For smaller A's, one obtains consecutive excited states of ([5]). This is confirmed 
by Fig. l(a,b). The taller and/or wider the potential, the more high-frequency modes with 
localized envelopes there exist. The non-localized envelopes, corresponding to AA < is 
([5]), still have spatial features on the scale of order one (see Fig. 2(a)). For an attractive 
potential in ([1]), V{x) < 0, there exist no localized solutions of ([5]). In that case, the 
envelopes of all high-frequency modes are not localized. Such envelopes of the first and fifth 
highest-frequency modes for V{x) = — 3sech(0.5x) are shown in Fig. 2(b,c). 




Figure 2: As in Fig. 1, solid and dashed lines show the moduli of the mode envelope 
and of the potential, (a) The fifth (non-localized) highest-frequency mode obtained by 
the finite-diff'erence approximation of ([T]) with V = 3sech(0.5x). (b,c) First and fourth 
highest-frequency modes (non- localized) for the potential V = — 3sech(0.5x). 
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